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Abstract: Scanning laser ophthalmoscopy has been used to measure 
individual cone-photoreceptor directionalities in the living human eye. The 
directionality is determined at different retinal eccentricities where it is 
expected that cones have diameters ranging between 5-10/im, comparable 
to the spot size of the incident beam. Individual cone directionality values 
are compared with the predicted directionalities obtained by using the 
waveguide model of light coupling to and from photoreceptors for the case 
of a focused incident beam. 
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1. Introduction 

The human eye is directionally sensitive as demonstrated by the Stiles-Crawford effect (SCE) 
[1] seen as a dramatic reduction in visibility of a light beam impinging the eye pupil towards its 
rim. The directionality also manifests itself in light reflected off the retina in what is referred to 
as the Optical Stiles-Crawford effect (OSCE) [2]. The SCE and OSCE are normally fitted to a 
Gaussian function at the pupil plane given by [1]: 

ij = io- p( - r - r ° )2 (i) 

where r is a point within the pupil, ro is peak position of the intensity distribution and p is the 
coefficient of directionality that measures the width of the pupil apodization. 

The SCE and the OSCE are related to the waveguide properties of the cones and their di- 
rectionality selective mechanisms [3], A well-accepted model of the ideal waveguide cone de- 
scribed by Snyder and Pask [4] assumes that the SCE is created by individual cone waveguide 
effects, where the light coupling and its consequences for the apodization effect at the pupil 
are collective contributions of individual cone effects. It has been experimentally observed that 
optical fiber properties of an ensemble of cones resemble those of an individual cone [5] and 
that the photoreceptors are sensitive to the slope of the wavefront incident on the retina [6] as 
expected on account of the SCE model. 

Vohnsen et al. [7] have upgraded the waveguide model introducing the diffraction properties 
of light to and from photoreceptors. Unlike a ray-optic geometrical approach, the diffraction 
model is needed to analyze the directionality and indeed an analytic p may be found including 
the cases of focused and non-focused illumination on the retina [7]. Actually, for a focused 
beam of light, the directionality turns out to be highly dependent on the ratio between the spot 
size of the incident light and the size of the mode propagated along the photoreceptor. The 
model has been examined experimentally using a microstructure fiber bundle to simulate the 
retina cone mosaic [8]. 

Different techniques have been implemented to measure directionality of cones using psy- 
chophysical and reflectance measurements. In the original Stiles-Crawford setup [1] and poste- 
rior variations as introduced by Enoch & Hope [9], the entrance beam angle is changed while 
the subject compares the reduction in visibility with that of a reference beam entering the eye at 
the pupil point of maximum visibility. Gorrand & Delori [ 10, 1 1 ] have developed a photorecep- 
tor alignment reflectometer that determines photoreceptor alignment by measuring the distribu- 
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tion of reflected light in the pupil. Delint et al. [12] introduced scanning laser ophthalmoscopy 
to measure local directionality in areas larger than 2 x2deg (i.e. s» 0.5 x 0.5mm). In both of these 
reflection techniques, two apertures are used to simultaneously alter the angle of the entrance 
and exit beam, scanning the subject's pupil. Burns et al. have modified the Gorrand and Delori 
setup avoiding to stop the exit beam at the pupil permitting an image of the entire distribution 
of light backscattered by the retina [13]. Gao et al. have explored optical coherence tomogra- 
phy to measure directionalities at different major layers of the retina [14], and posteriorly have 
implemented the Hartmann-Shack wavefront sensor to distinguish the guided (directional) and 
the diffused (non-directional) components of the reflected light from the retina [15]. Exploring 
the gain of resolution that adaptive optics provides, Roorda & Williams have measured indi- 
vidual human cone directionalities using a fundus camera for a retinal region 1° nasal from the 
fovea [5]. 

Here, individual-human cone directionality measurements carried out with a scanning laser 
ophthalmoscope (SLO) are presented. The directionality coefficients are obtained at different 
superior parafoveal eccentricities of the authors's right eyes where the size of the cones are 
expected to increase from approximately 2.5/im to ^10/im at 1° to 20° from the fovea, respec- 
tively [16]. The values are analyzed with the model of light coupling to and from photoreceptors 
considering an incident spot size of w r w 5.2/im (i.e., a 10.4/im beam waist). The individual 
directionalities are compared with the corresponding collective values and directionality values 
obtained in previous studies. 

2. Experimental method 

Retinal images are acquired at a frame rate of 47fps (frame resolution of 512x490 pixels) using 
a non-commercial SLO. The SLO contains vertical and horizontal scanning mirrors located at 
conjugated planes of the subject's pupil (P), a (2x) beam expander, and an adaptive optics 
system (reducing the total optical aberrations to a maximum remnant RMS wavefront error of 
approximately X /4). Experimental details about the SLO system are contained in Ref. [17]. 
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Fig. 1. Experimental method, (a) Experimental setup and (b) CCD images of the beam 
profile at the conjugated (P') and subject's (P) pupil position. * Further details about the 
SLO are given in Ref. [17]. 



As represented in Fig. 1, a laser diode at 785 nm wavelength (and 1 nm bandwidth) has 
been used as source. A laser diode was chosen for convenience as it was already present in our 
laboratory and in use for a related study [17], although to reduce unwanted interference effects 
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a superluminescent diode would be preferable. Nonetheless, the images and videos obtained 
(Section 4) do confirmatively show that cones have been resolved. An aperture I\ with diameter 
of ~ 0.8mm is translated 3.5mm through the conjugated pupil's plane (P') by using a motorized 
stage, corresponding to a movable artificial pupil ) of ~ 1.6mm translated 7mm across the 
plane of the subject's pupil. A beam splitter (BS) is used to select a part of the reflected light 
which is filtered by a confocal pinhole (C-P) and detected by an avalanche photodiode (APD). 
100 image frames are collected simultaneously with the motorized transversal of 7mm of the 
incident illumination. Thus, between each image frame the incident beam has been displaced 
70/im corresponding to a full scan of the subject's pupil in a lapse time of 2 seconds. 

As shown in Fig. 1(b) the aperture scan through the laser introduces a Gaussian intensity- 
profile, where 0.9 mW (r=0) and 0.3mW (r=+3.5mm) are the maximum and minimum incident 
powers at the subject's pupil. The Gaussian profile of the laser has been characterized at the 
pupil plane in steps of 0.1mm and the Gaussian fitting-curve parameters are used to correct the 
reflection intensity in the videos such that the illumination corresponds to an intensity-constant 
artificial pupil. Figure 4(b) (Media 3) shows an example of an unprocessed and processed pair 
of videos, i.e., without and with correction of the beam intensity variations. 

The right eye of both authors (DR: age 30 and BV:41, both with normal vision) has been 
examined with the dilated pupil (« 8mm size, applying tropicamide 1%). A fixation target 
(Green LED) was used to set the viewing angle respectively at 2.5°, 5°, 10°, 15° and 20° 
superior thereby probing the superior parafoveal retina. Ethical approval for the research was 
obtained from the UCD Human Research Ethics Committee. 

The directionality values are obtained acquiring the intensity values within the selected area 
(using ImageJ software) in steps of 0.35mm (i.e. every 5 frames) through the pupil scan (corre- 
sponding to 20 experimental points), obtaining an accurate fit of the apodization effect by using 
Eq. (1). Under unparalyzed accommodation conditions and during a time-lapse of 2 seconds 
there are maximum RMS wave-front error variations of a ^0.02jUm with random character- 
istics [18], corresponding to a random Strehl ratio variations of A5=0.02. Thus, it is expected 
that in our case with a paralyzed accommodation condition, wave-front variations cannot alter 
significantly the Gaussian apodization of the reflection intensity during the scan of the subject's 
pupil. 

3. Waveguide model of light coupling [7] 

As represented in Fig. 2, the light guided in the cone is propagated from the inner to the outer 
segments where it is reflected by backscattering processes. Posteriorly, the light is guided back 
predominantly axially with single-mode characteristics (Gaussian shape), diffracted at the exit 
of the photoreceptor, and creating a Gaussian intensity distribution at the subject's pupil. For 
the case of an SLO system a focused beam is used to scan the subject's retina, in our experiment 
the estimated spot size at the retina is w r « 5.2/im for an aperture of I\ ~ 1.6mm at the pupil. 
The incident beam intersects the photoreceptor at an angle 6 = rj ' f eye where part of the light 
is coupled to discrete modes as governed by waveguide boundary conditions. The fraction of 
power transmitted to the photoreceptor if incident with its peak value at the photoreceptor axis 
can be estimated from [19]: 

-2(CT eve w r w„,) 2 0 2 "| 

where w m is the width of the mode propagated. The eye parameters used are those of the 
schematic eye, i.e., a focal length / e>f ,=22.2mm and a constant index of refraction n f5 , f .=1.33. 

Assuming that the light intensity reflected by the photoreceptors towards the pupil [Eq. (1)] 
is linearly related to the amount of coupled incident light (i.e. bleached or constant fraction of 
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absorption) a directionality coefficient can be inferred directly from Eq. (2). as: 

2 
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Fig. 2. Schematic of light coupling to the cone-photoreceptors at the retina plane for an 
SLO system. 

Because of the small size of the outer segment and the scattering mechanisms within the cone 
it is expected that only the fundamental mode of the waveguide is guided back [7]. For a single- 
mode condition the mode-width is related to the radius of the waveguide (i.e. the photoreceptor 
radius, a) given by the empirical relation [20]: 
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where V is the waveguide parameter. In Fig. 3 the dependence of the directionality coefficient 
with the photoreceptor size is estimated for different incident spot-size conditions and a waveg- 
uide parameter of V=2.Q by use of Eq. (3) and Eq. (4) [8]. 
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Photoreceptor radius, a (urn) 

Fig. 3. Dependence of the predicted directionality coefficient with the photoreceptor size. 
The values are estimated for an illumination spot-size of w,-= 2.5/im, 3.5/im and 5.5/im . 



Due to the scanning, only an area given by the spot size of the illumination beam (in our 
case approximately 1-3 cones) is illuminated simultaneously at any given instant and possible 
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collective effects as scattering created by a distribution of cones is avoided. This is under the 
assumption that multiple scattering may be ignored, which appears reasonable considering the 
low contrast of refractive indices within the retina whereby the illumination is dominated by 
the focused incident beam. 

4. Results and discussion 

In virtue of the small size of the cones in the central region, an individual cone analysis is not 
well-controlled, limited by the resolution of the image, and instead average values of selected 
areas are determined. Figure 4(a) (Media 1) and Fig. 5(a) (Media 6) show the directionality 
values obtained in areas selected at a distance of 0.65°, 1° and 1.5° from the center of the fovea 
for both subjects. 




Fig. 4. Retinal images of the subject BV and directionality values at different retinal re- 
gions: (a) eccentricity from 0° to 1.5° (Media 1), the blue zones are the areas selected 
to measure the averaged value of directionality, the directionality coefficient of the areas 
selected is given in mm -2 units; (b) superior retinal eccentricities of 5° (Media 2), 10°* 
(Media 3), 15° (Media 4), and 20° (Media 5). The individual directionality values are rep- 
resented by: (Yellow box) p= 0.10 - 0.15mm" 2 , (Red box) p= 0.15 - 0.20mm" 2 , (Yel- 
low star) p= 0.20 - 0.30mm" 2 , (Red star) p > 0.30mm" 2 . The opened square represents 
the area selected to calculate the directionality averaged values with (Yellow square) p= 
0.10 - 0.15mm" 2 and (Red square) p= 0.15 - 0.20mm" 2 . * The media shows a processed 
(intensity-constant artificial pupil) (a) and an unprocessed (uncorrected for the Gaussian 
laser beam profile) video (b). 

Taking into account the radial eccentricity dependence of the cone-size [16], a ring-shape for 
the area selection has been adopted to analyze the directionalities. The highlighted blue regions 
are the areas selected for the analysis, a histogram furnish the mean intensity of the pixels in 
the area for each frame using an ImageJ application. Software-tracking and averaging of frames 
have been avoided, providing non-altered intensity values of reflection. Near the fovea there is 
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a brighter zone that follows the scan of the pupil caused by reflection at the inner limiting 
membrane [21]. It has a little impact, however, on the Gaussian shape of the apodization effect. 
The same effect cannot be seen for the parafoveal regions due to the changed direction of gaze. 

As shown in Fig. 4(a) and Fig. 5(a) the increase of the directionality with the eccentricity 
resembles the dependence reported by Burns et al. [22] for OSCE directionality measurements, 
and the values reported by Westheimer [23] and Enoch & Hope [9] for the SCE measurements. 




Fig. 5. Retinal images of the subject DR and directionality values at different retinal re- 
gions: (a) eccentricity from 0° to 1.5° (Media 6), the blue zones are the areas selected 
to measure the averaged value of directionality, the directionality coefficient of the areas 
selected is given in mm -2 units; (b) superior retinal eccentricities of 5° (Media 7), 10° 
(Media 8), 15° (Media 9), and 20° (Media 10). The individual directionality values are rep- 
resented by: (Yellow box) p= 0.10 - 0.15mm -2 , (Red box) p= 0.15 - 0.20mm -2 , (Yellow 
star) p= 0.20 - 0.30mm -2 , (Red star) p > 0.30mm -2 . The opened square represents the 
area selected to calculate the directionality averaged values since (Yellow square) p= 0.10 
- 0.15mra -2 and (Red square) p= 0.15 - 0.20mm -2 

Individual directionality measurements are possible to obtain with the method at eccen- 
tricities larger than 2°. Figure 4(b) (Media 2-Media 5) and Fig. 5(b) (Media 7-Media 10) 
show images and videos taken at different retinal eccentricities (5°, 10°, 15° and 20° superior 
parafovea). The positions of a series of cones have been selected randomly across the cone 
mosaic and the average reflected intensity for each cone has been measured integrating across 
a circular region with diameter corresponding to 8 pixels (i.e., approximately 4 micron) in each 
frame. Small movements of the subject's eye change subtle the image location of the cones so 
that in each case the position of the selected area is shifted manually (frame by frame) tracking 
the center of the cone analyzed. As shown in Fig. 4(b) and Fig. 5(b) the symbols (■ and *) 
represent the individual directionalities. The opened squares correspond to the directionality 
value of the area enclosed by each square. 

The individual directionality values range from 0.08 to 0.12mm -2 at 2.5° of eccentricity and 
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increase approximately three times at 20° where it is expected that cone diameters are similar 
to the spot size of the incident beam. However, if the directionality is measured for areas across 
a group of cones the directionality is about 0.14 mrrC 2 (represented by opened squares in Fig. 
4(b) and Fig. 5(b)). If the total area of the image is analyzed the directionality is approximately 
0.12mm -2 , and approximately constant for different eccentricities. 

Figure 6(a) shows the averaged apodization effect of the full images for the subject (BV) for 
the retinal eccentricities analyzed, where the peak position and the directionality are approxi- 
mately constant. The peak positions of the apodization effect are approximately equal for both 
subjects (for simplicity not shown here) and have been reported previously for psychophysical 
Stiles-Crawford measurements [25]. 
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Fig. 6. Reflected average image intensity as a function of pupil entrance point of the imag- 
ing beam for the subject BV: (a) total average for superior retinal eccentricities of 5°, 10°, 
15°, and 20°, the error values represent the standard deviation; (b) Two examples of in- 
dividual cone values compared with the value obtained analyzing an area containing the 
same cones for an eccentricity of 10°, the error values have not been represented such that 
variations in peak locations can be easily observed. 



The difference between the averaged and individual directionality values may plausible be 
explained by the following reasons: the existence of rods occupying the inter-cone spaces for 
eccentricies larger than 3 degrees [16], where directionality values are smaller than those of 
cones [24], and/or a difference in pointing direction for each cone that when averaged creates 
a broader directionality effect [5]. As shown in Fig. 6(b), the peak position of the apodization 
effect for individual cones shows a subtle difference of pointing direction that can introduce a 
broader directionality effect when the area containing the same cones is analyzed. 

The directionality values obtained for both subjects at the studied retinal regions are shown 
in Fig. 7. For the total average values, all the image area is analyzed, and the directionality is 
obtained for a total of 7 videos. Using the retinal mosaic images at different retinal eccentricities 
reported by Curcio et al. for ex-vivo human eyes [16] (i.e. from the fovea up to an eccentricity 
region of 16mm), it is possible to estimate cone diameters at a given eccentricity. A polynomial 
fitting is used to describe this dependence that by means of Eq. (3) results in the solid line 
shown in Fig. 7. 

In the particular case of Fig. 5(b) (Media 7) the reflection of the blood vessel appear with 
similar brightness as the cones, indication a non ideal focus on the photoreceptors. For the 
subject DR it suggests an imperfect conjugation with the confocal pinhole or an incomplete 
defocus correction. Despite of this, defocus variations are smaller than AZ5=0.05/im that would 
increase the spot size at the retina plane from 5.2/im to « 6-7 /im, hence the experimental 
directionality dependence with eccentricity (solid line in Fig. 7) is still well predicted by use 
of Eq. (3). Furthermore, it can be noted that the blood vessel reflection has not a directional 
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Fig. 7. Individual and total area-averaged directionalities of both subjects at different supe- 
rior parafoveal eccentricities. The solid line corresponds to the predicted dependence of the 
directionality coefficient with the eccentricity using Eq. (3) for the special case of a = w m 
using the retinal mosaic reported by Curcio et al. [16]. 



characteristic appearing approximately constant during the pupil scan. 

As can be seen in Fig. 5(b) (Media 8, Media 9), there are variations of reflected intensity un- 
correlated with the intensity variations of neighboring regions, usually present for a few frames 
(1-10 frames). The flickering effects can be easier identified for videos at retinal eccentricities 
larger than 10° and are repeatable for a trial of videos for the same retinal localization. How- 
ever, most regions with light flash characteristics are smaller than the spatial resolution of our 
SLO, and the intensity variations have not a well-defined apodization shape. The origin of the 
effect is not fully understood at present but may be explained for the following reasons: the 
coherent light illumination is causing the changes either due to rapid index/length variations or 
local intensity speckles that modulate the coupling efficiency, and/or some cones may have a 
very narrow acceptance angle because of their composition and geometry. 

5. Conclusion 

We have demonstrated an SLO imaging method to measure individual cone-photoreceptor di- 
rectionalities. The method explores the resolution of the SLO system which combined with 
beam-scan at the pupil allows for a fast and accurate acquisition of the photoreceptor accep- 
tance angle. 

As demonstrated by theory and experimental values, the directionality is highly dependent 
on the light coupling mechanism of the incident beam to the photoreceptor, in this case partic- 
ularly manifested in the relation between the spot size of the incident beam and the diameter of 
the photoreceptor. We have observed that although the coupling of the individual cones to the 
focused incident beam is highly directional a smaller directionality is found when averaged on 
a larger retinal area. 
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